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SUMMARY 

The degree of approximation used in determining the 
orbits of earth satellites is reflected in the residuals (dif- 
ferences between calculated and observed positions). The 
least- squares procedure generally used to fit theory to ob- 
servation tends to obscure the significance of theoretical 
parameters, so that the physical sources of residuals cease 
to be apparent. A method is outlined herein for estimating 
the magnitude of the residuals to be expected from an ap- 
proximate theory presumed to have one missing or incor- 
rect term. 
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A METHOD OF ESTIMATING RESIDUALS 
IN ORBITAL THEORY 


INTRODUCTION 

fe determining the orbits o, earth satellites, as 
approximate theory is used. The degree o Analysis oi «te re- 

££ “ ^ “esent understanding o, the pertinent physical 

phenomena. 

This paper outlines a method 

caused by certain inadequacies in p ^ least-squares minimization proce- 

ci:rr,:r:rrr; : T 

"r:r:r“"^^^ate «. ... - 

provide a model for the calculations to follow. 

Ior : re^r^^ * sivcn 

rigorously by 

( 1 ) 


s = 


S 0 + Vt + 


vt- 


( 2 ) 


• * rxf +hP v term had not been known, and that the ob- 
Suppose further that the existence of the V term naa 

ser vat ions had been fitted to 

r :”sup"ts is ^ 

r , it might first be hypothesized that R is of the order 
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R " T / (S " S') 2 dt = AS 2 + + As Av T + ASVT 2 AVVT 3 

o j zv 34 

where 

AS = S o ■ s o • Av = V - V' . 

In the symmetrized form, let 

V = x 0 , AS = X j , and AV - x 2 . 

Then 
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Note that in tins symmetrized formulation, may always be written as the ditference 
of two terms. The first term represents the residuals introduced by die missing or in- 
correct terms in the approximate theory (Equation 3). The second term represents the 

;::r t?rf ,alS e “ eCte<i by lhe mlnimlZa “°" I™—- This term measures the 
e«ent to which the free parameters can compensate for the incorrect or missing param- 
eter (and in the process separate themselves from their direct physical analogs). 
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ERRORS DUE TO DRAG 

The primary atmospheric effect (i.e., drag) is a secular shortening of the orbital 
period. For the satellite 1958 Beta 2, Jacchia* found that the time derivative of the 
orbital period was sinusoidal with time, with a period of about one month. This suggests 
that a constant second time derivative of the orbital periods may be used for arcs less 
than one week, while for longer arcs the sinusoidal variation should be represented. 

In this section it is assumed that the approximate theory allows the period to vary 
linearly with time; and the residuals due to neglecting these terms are estimated. 

From the discussion on observables and the minimization of residuals, it is seen that 
R Cot') 2 + (i-i) 2 ; [oM 2 cos 2 0 + cm 2 ] ( 10 ) 

If the assumptions are restricted to circular orbits in a central gravitational field, then 

R = jj (Sf) 2 , (11) 


where f is the true anomaly. 

Two calculations are now made: 
a) n is quadratic ally time dependent : 
The "true" f is given by 


f f n ^ nt + 4- nt 2 4 i- nt 3 ; 

u 2 6 

(12) 

while the computed f is taken as 


f ' ^ n't + y n't 2 , 

(12a) 

where n is the mean motion. The minimization of R with respect to f^, n', 
the condition that 

, n ’ exacts 

+ 120 ’ 

(13) 

n ' - n ~ ^ nt 2 , 

(13a) 

n ' = n + ~ iit , 

(13b) 


^Jacchia, L., Nature 183(4660): 526-527, February 21, 1959; and 183(4676): 1662 - 1663, June 
13, 1959 
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h<t> = x 0 


jnt (X - 0 2 ) (4^ 2 “ 1) cos 0 
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With the foregoing, construct 

t 

r = -i-V a ij x i x j = y f ( s ^) 2 ( Sm2 ) dt 

P i j •'n 


(17) 


(18) 


(19) 


If the time average over an integral number k of orbital periods is taken, it is found that 


a„„ = Kr N 2 0, + 0(j 2 ) 
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and 


a oi " a io 


jN 2 g 3 + I^oo 

3 3 


a 02 3 20 


j 1 2 n 2 0 5 

+ 0(j) , 


a 0 3 a 30 
a 04 ~ a 40 

a 1 2 = a 21 
a 1 3 = 3 31 

3 1 4 = a 4 1 

a 23 ~ 3 32 

a 24 ” a 42 

3 34 = a 43 

where 
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In general, it seems necessary to solve for R min numerically. However, there are 
two special cases of interest that are tractable algebraically: 

(1) Setting Xj = x 2 = 0 and thus minimizing only with respect to the initial values of 

the true anomaly and the longitude of the node, it is found that 


= 


*4 = 
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2 • 2 « T 2 , , 2 
x 0 J N 

12 


( 20 ) 


This result is consistent with that usually found when compensating for an error term 
linear with time by correcting the initial position. 


(2) Setting x 3 = x 4 = 0 and letting k > 100 , then 


where 


X 0 2 j 2 N 2 2 
3 



*3 2 + i 2 w 2 ^g x \ 
1 + j 2 N 2 ^2 J 


e 


2 



( 21 ) 


In general, the correction term here is quite small; this indicates that the mean motion 
and the inclination cannot compensate for an incorrect j. 


CONCLUSION 

In this paper a general procedure has been outlined for estimating the magnitude of 
the residuals expected from an approximate theory presumed to have one missing or in- 
correct term. If, as is usually the case, a least- squares procedure is used to fit theory 
to observation, a reasonable estimate of the residuals cannot be obtained by considering 
only the missing terms. Rather, the remaining free variables may compensate for the 
missing terms, to reduce the residuals by an order of magnitude. 
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